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Abstract 

In this article, we study the self-similar solutions of the 2-component Camassa-Holm equa- 
tions 

pt + upx + pux = 

rat + 2uxm + unix + appx = 



(1) 



with 



m = u — a Ux 



(2) 



By the separation method, we can obtain a class of blowup or global solutions for a = 1 or 
— 1. In particular, for the integrable system with cr = 1, we have the global solutions: 



p{t,x) = < 



a{s) 



fM for 71^ < -^ 



l(t,X) — ^.gj. 



0, for V >Y- 
^ - 0, a(0) = ao > 0, a(0) = ai 



(3) 



3a(s)l/3 



where n = , ^-, ,„ with s = 3t; f > and a > are arbitrary constants. 
Our analytical solutions could provide concrete examples for testing the validation and stabil- 
ities of numerical methods for the systems. 
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1 Introduction 

The 2-coniponcnt Camassa-Holm equations can be expressed in the following form 

pt + upx + pux = 

(4) 
mt + 2uxm + unix + crppx — 

V 

with 

m = u- a^Uxx, (5) 

X G R and u = u{x, t) E R is the velocity of fluid and p ^ p{t, x) > is the density of fluid. The 
constant a is equal to 1 or —1. If cr = —1, the gravity acceleration points upwards [2], [3], |13j . 
[IB] and [TS]. If cr = 1, there exist some papers regarding the corresponding models [5], [H], [IB] 
and [13]. When p = 0, the system returns to the Camassa-Holm equation [T]. The searching of a 
model equation which can capture breaking waves and peaked traveling waves is a long-standing 
open problem |24j . Here, the Camassa-Holm equation satisfies the above conditions as a model 
equation. The searching of peaked traveling waves is motivated by the wish to discover waves 
replicating a characteristic for the wave of great height (waves of largest amplitude), that are exact 
traveling solutions of the shallow water equations, whether periodic or solitary [S], [23] and [5]. 
The breaking waves can be understood by solutions which remain bounded but the slope at some 
point becomes unbounded in a finite time [6]. Meanwhile, there is an alternative derivation of the 
Camassa-Holm equation in [TH] and [TU] . 

With a ~ I, the system ^ is integrable [T7] and [T]. It can be expressed as a compatibility 
condition of two linear systems (Lax pair): 

V'.. = (C'p' + Cm + J)V (6) 
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and 



■0* = (^ - u)tpa: + -U^ll) 



with a spectral parameter C,. 

It is bi-Hamiltonian, the first Poisson bracket 



{^1,^2} = - 



- — imd + am)- h - — pd— h -^op- — 

dm dm dm op dp dm 



with the Hamiltonian 



and the second Poisson bracket 



dx 



H = — (mu + p ) dx 



(7) 



(8) 



(9) 



{Fi,F2} = -J 



5Fi ^dF2 , 5F.dF2 

dm dm dp dp 



dx 



with the Hamiltonian 



There are two Casimirs: 



and 



H = - {up^ + M^ + uul) dx. 



Mass = / pdx 



m,dx. 



(10) 



(11) 



(12) 



(13) 



In this article, we adopt an alternative approach (method of separation) to study some self- 
similar solutions of 2-component Camassa-Holm equations Q- Indeed, we observe that the isen- 
tropic Euler, Euler-Poisson, Navier-Stokes and Navier-Stokes-Poisson systems are written by: 



Pt+V ■ (pu) =0 
p[Uf. + (u-V) u)] + VP = -pV^ + vis{p, u) 
A<^{t,x) ^a{N)p 

where a{N) is a constant related to the unit ball in R^ : a{l) — 2; a(2) = 2tt and for A^ > 3, 



(14) 



a{N) = N{N - 2)Vol{N) = N{N - 2) 



r(iV/2 + 1) 



(15) 
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where Vol{N) is the volume of the unit bah in R^ and F is a Gamma function. And as usual, 
p = p{t, x) and u = ■u{t, x) £ R^ are the density and velocity respectively. P ~ P{p) = Kp^ is the 
pressure, the constant K >Q and 7 > 1. And vis{p,u) is the viscosity function. 
We may seek the radial solutions 

p{t, x) = p{t, r),u^ -V{t, r) =: -V , (16) 

r r 

with r ^ (j2Zi ^1) ■ 

By the standard computation, the Euler equations in radial symmetry can be written in the 

following form: 

Pt + Vpr + pVr H pV = 

(17) 
For the mass equation in radial symmetry, (jl7P i . we well know the solutions' structure (Lemma 3, 

my- 



f'^W)) ,, . a{r) 

— , , ., , u(t,r) = , ', 

a(i)^ a{r) 



P(t'^)--J7Sr^ "(*'0=77l(^- (18) 



As the 2-component Camassa-Holm equations, ((4]), are very similar to the Euler system p7)) . 
in some senses, we can apply the separation method ([2], [21], [20], [25], [27]) to the systems 
(UJ. In fact, we can deduce the nonlinear partial differential equations into much simpler ordinary 
differential equations. In this way, we can contribute a new class of self-similar solutions in the 
following theorem: 

Theorem 1 We define the function a{s) is the solution of the Emden equation: 

a(0) = flo 7^ 0, a(0) == a^ 
and 



(19) 



«")^!fi"'+(?)^ ,20, 

where rj — , f^^a with s — 3t; a > 0, ^ ^ 0, Oq and ai are arbitrary constants. 

For the 2-compenent Camassa-Holm equations Bl), there exists a family of solutions: 

(1) for a = -1, 
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(la) with ^ < and ao > 0, 



^(''^) = {^^' ^} '"(''") = ^^- ('') 



The solution i21\) blows up in a finite time T. 
(lb) with ^ > and oq < 0, 



"('-' = il^' "('-> = §|- (^^' 



The solution i22\) exists globally; 

(2) for a = 1, 

(2a) with ^ > and ao > 0, 



.(...). ma. {-^,0},.(,,x).g|x^ ,23) 



The solution h23\) exists globally. 
(2b) with C < and oq < 0, 



^<'-' = il^' "('-> = SI- (^*> 



The solution {24-^ blows up in a finite time T . 



2 Separation Method 

First, we can design a nice functional structure for the mass equation: 
Lemma 2 For the 1-dimensional equation of mass OPi.' 

Pt + upx + pu^, = 0, (25) 



there exist solutions, 



,. . fL(3ty/-^^ ,, . a(3t) 



with the form f {ri) > E C^ with rj = 'tSyTJI j '^'^'^ a(3t) > G C^ 
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Proof. We just plug ^ into ^ to check: 

Pt + upx + pux 

_ d M (a(3t)i/3) \ a(3t) d f fia{3ty/3)\ /(a(3t)i/3) 9 / a(3i) 
~ 9f 1^ a(3t)i/3 j ^^py'^a^ [ a(3f)i/3 j a(3i)i/3 ft^ ^^pj' 

a 
(3?) 



1/3+1 ( 3) • «(3i) • 3 • fi^r^^^l/3)+ ai'StW^^^ a(3tW^^ dtLiitV/^^ 



a(3t)i/3+i^ 3' ' ' •''a(3t)i/3^ a(3t)i/3'' ^a(3t) 
a(3i)a; /( 0(34)1/3) d f x \ fiaisty/^' a{3t) 



(3t) a(3t) ax Va(3t)i/V a(3t)i/3 a{3t) 



-1I/3+1 ^ „("^/U/3-' S 



WTQ°(3i)3 



a(3t)i/3+i a(3i)i/3 a(3f)i/3'a(3t)i/3+i3 

a(3^)a;/(a(3f)i/3) 1 /(a(3t)i/3) a(3t) 

a(3t) a(3i)i/3 a(3t)i/3 + a(3t)i/3 a{3t) 

= 0. 

The proof is completed. ■ 

In [11] and [26] , the qualitative properties of the Emden equation 

a(0) = ao > 0, a(0) = ai, 
v 

where iV > 2, 

were studied. Therefore, the similar local existence of the Emden equation HHJ), 



(27) 
(28) 

(29) 

(30) 

(31) 

(32) 
(33) 



(34) 



(35) 



a(0) = oo 7^ 0, a(0) — ai, 

can be proved by the standard fixed point theorem [11] and [26] . To additionally show the blowup 
property of the time function a{s), the following lemmas are needed. 

Lemma 3 For the Emden equation I119\) . 



a(0) = flo > 0, d(0) = ai, 



(36) 
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(1) if ^ <0, there exists a finite time S, such that lim a{s) = 0. 

(2) if ^> 0, the solution a{s) globally exists, such that 



lim a{s) = +00. (37) 

s— >-+oo 

Proof. (1) For the Emden equation ff36|l . we can multiply d(s) and then integrate it, as follows: 

(38) 



^ - l^aisf/^ 



with the constant 

for C < 0. 

The equation (|38| becomes 



Oi 3 2/3 ^ f, fr,r,\ 



d(s) = ±^2e + 3^a{s)^/^. (40) 

We can define the kinetic energy thus: 

Fk^n ■■= ^, (41) 



and the potential energy thus: 



The total energy is conserved thus: 



Fpot = -l^aisf/'. (42) 



-r{Fk^n + Fpot) ^ 0. (43) 

ds 
By the classical energy method for conservative systems (in Section 4.3 of yS]); the solution has a 
trajectory. The time for traveling the orbit [0,asup] can be estimated: 

S^rds<2 r ""^'^ - 2 /^^'" '^^^^ (44) 

^0 Ja..,, ^2(0+M£(£)!^) Jo ^2(0+3SaM!^) 



Then we let 6(5) := J^a(s)i/3^ j^^ve 



dGis) = /J^^, (45) 

V 6 a(s)''''^ 
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and 



da{s) 



a{sf^^dGis) 



-3CV C 



G{sfdG{s). 



(46) 



Then, the equation (HH becomes: 



f2B\3/2 
y 3() 



da{s) 



2V2 P6 



Ve 



G{s)^dG{s) 



^/o(g I 3Ca(s)V3 - -SeV e 7o v/(^-G(s)2)' 

The integral is solvable exactly by the integration formula 17.11.3 on page 82 of [22] : 



(47) 



^ x'^dx 



V(^^^ 



?;V^ 



xy/tf — X 



■ -i^\ 

- sm —= 

2 vv^y 



x=\/e 






(48) 



Therefore, we showed that there exists a finite time S, such that lim a{s) — 0. 

s->S- 

(2) for the case of ^ > 0, the potential function is 



F, 



pot 



-^^ais)^'- 



As the infimum of a{s) is: 



,-261,3 
Oinf = (^)^ >0, 



(49) 



(50) 



the solution a{s) is uniformly bounded below. The time for achieving the infimum is finite if ai < 0: 

"^1 '■"'' da{s) 



Si^ 



Jo J(^) 



< +00. 



(51) 



Therefore, for any constant Oi, the solution a{s) must increase, after some finite time. And the 
time for traveling the interval [inf a(s), +00) or [ai,+oo), is infinite: 

da{s) 



So^ 



+ 00 



(^f)' a/2(0 + 



3ga(a)^/3 



= +00. 



(52) 



Therefore, we showed that the solution a{s) globally exists, such that 

lim a{s) = +00. 

The proof is completed. ■ 

We may let b{t) — ~a{t), in the Emden equation (|19p . for oq < 0: 

a(0) = flo < 0, a(0) = Oi, 



(53) 



(54) 
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have 



ks) 



J_ 







3h(s)i/3 

b{0) = -ao > 0, 6(0) = -ai. 
Alternatively, the ordinary differential equation can be rewritten as the dynamical system: 



(55) 



£^is) = vis) 
Ay(s) - « (56) 

a;(0)-ao^O, y(0) = ai. 

It is symmetric about the origin (0,0). Therefore, we can use the above lemma to drive the 
corresponding lemma for oq < 0: 



Lemma 4 For the Emden equation ^_ 



a(0) = ao < 0, d(0) = ai, 



(1) i/^ < 0, there exists a finite time S, such that lim a{s) — 0. 

(2) if £, > 0, the solution a{t) globally exists, such that 



lim a{s) — — cx). 

s— ^+oo 



(57) 



(58) 



After obtaining the nice structure of solutions (|26p . we simply use the techniques of separation 
of variables ([W], [21], [20], [25] and [27]), to prove the theorem: 

Proof of Theorem [1] From Lemma [2 it is clear to see our functions ([21]), ([22]), ^^ and (p4)) . 
fit well into the mass equation, (|3])i, except for two boundary points. 

The second equation of 2-component Camassa-Holm equations (01)2, becomes: 



rut + 2uxm + urrix + <JpPx 
= (u - Uxx)t + 2u2:(u - Uxx) + u{u ~ Uxx)x + crppx- 



(59) 
(60) 



As our velocity u, in the solutions (HH, ([2^ . (P5|) and (IMl) . is a linear flow: 



a(3i) 

u = ,„ .s X, 



a(3t)' 



(61) 
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we have 



The equation (|60| becomes 



u^x = 0. (62) 



ut + iuxu + appx (63) 

d fd{3t)\ „/a(3i)\ a(3t) , fiaiaty/^) ( fiaisty/^)] 



9t \a{3t)J \a{3t)J a{3t) a(3t)i/3 I a(3i)i/3 

,a(3t) X3t)2 N^ /a{3t)\a{3t) /(^Tisfyv^) /(^(#7^) 1 ,„_. 



a(3t) a(3i)V Va(3t)/ a{3t) a(3t)i/3 a(3t)i/3 a(3t)V3 

a(3<) •/ (a(3t)i/3)/(a(3t)i/3) 

= 3— — x + cr ^^ — (66) 

a(6t) a[6t) 

"" '^^V + fMHv)), (67) 



«(«) - ?wlu7^ = 



a(3i) \cr 
with the Emden equation: 

" " - 3a(5)i/3 

a(0) = oq 7^ 0, d(0) — fli, 
by defining the variables s := 3i and rj := a;/a(s)^''^. 

Now, we can separate the partial differential equations into two ordinary differential equations 
Then, we only need to solve for - < 0, 



/(0) = -a<0; 
or for T > 0, 

/(0)=«>0; 
The ordinary differential equations ([M)) or (1701) can be solved exactly as 



(69) 



(70) 



/(.)^^Jf.^.(f)^. (71) 



In fact, we have the self-similar solutions in details: 
(1) for a = -1, 
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(la) with C < and oq > 0: 

p{t,x) = < 



/('') for ri^ < -s! 



_ a(3t). 



u{t,x)^ -^X 



a(3t)i/ 

0, for r/2 > -i|: 
"(*) ~ 3Q(f)i/3 = 0, a(0) =: ao > 0, a(0) = ai 



(lb) with C > and ao < 0: 



^(*'^) = ^W7^' «(^'2;) = - 



o(3t), 
0(3*7 '' 



°^(^) ~ 3a(f)i/3 = 0, a(0) = ao < 0, a(0) = ai 



/W=-eJt^2 + (|«)^ 



(2) for cr= 1, 

(2a) with C > and oq > 0: 



p(t,x) 



a(3t)i/3 



for 7/2 < ^ 



u(i, x) 



0, for if>^ 



a(3t) ^ 
a(3t)' 



(72) 



(73) 



(74) 



°^(*) ~ 3a(f)i/3 = 0' a(0) = ao > 0, a(0) = a^ 
(2b) with ^ < and Oq < 0: 

«(^) - M^ = 0' «(0) = ao < 0, a(0) - ai (75) 

With the assistance of Lemmas [3] and SI the blowup or global existence of the solutions can be 
determined under the prescribed conditions in the theorem. 
The proof is completed. ■ 



Remark 5 For ^ < 0, the blowup solutions i21\) and \24^ collapse at the origin 

lim pit, 0) = +00, 

with a finite time T; 

For ^ > 0, the global behavior of the solutions i22\) and Ii23\) at the origin is: 



(76) 



Km p{t, 0) = 0. 

i— )-+oo 



(77) 
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Remark 6 The solutions \21]) and \2S\) are only C^ Junctions, as the function /(r;) is discontin- 
uous at the two boundary points, for a > 0; 

lim /(77) ^ 0. (78) 

Remark 7 Our analytical solutions could provide concrete examples for testing the validation and 
stabilities of numerical methods for the systems. Additionally, our special solutions can shed some 
light on understanding of evolutionary pattern of the systems. 

Remark 8 For the integrable system with a ^ I, we may calculate the mass of 
(1) the blowup solution [24\ ) (or i75\)), £, < and ao < 0; 



Mass = ' "'" -^''- 



piO,x)dx = ^j^ / \ -7 — 
(2) the global solution ifSgj) (or (74\)), S, > and ao > 0; 



1 / a; \ /I ^ ^ 



a ) dx — +CX); (79) 



r+oo 
Mass^ I p{0,x)dx (80) 



-;''"i/| ^ ^ 



2 / \3 dx (81) 

"'o c ' 

2^ /-a?"/! 



1/3 



1 / X \ /a 



\ i \^J " VI 



2 

dx (82) 



2^/ y-^-^'+(^^) ds (83) 



« / 1 o ( Oi^ ^ 



2£ /""Vs a2 
"^ ' ' s^ds (84) 



V?7o V ^ 



2 

v!' 



(85) 



3 Blowup Rates 

We are interested in how fast the blowup solutions tend to the infinite as the time tends to the 
critical time T. The blowup rate of the constructed solutions (PT|) and ([M)) at the origin, is 
estimated in the following theorem: 
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Theorem 9 The blowup rate of the solutions i21\) and \24^ for a > 0, is 

lim_p{s,0){S-s)^/^>O{l). (86) 

s—^S 

Proof. We only need to study the blowup rate of the Emden equation (ITO|) , 



a(0) = ao > 0, d(0) = oi 
with ^ < 0. 
For the total energy 6 > 0, with the assistance of the equation pO)) . 



(87) 



S = s+ dr] = s+ -y-da ^ s- dn. > s + ' ^" 



la{s) da ^(,) ^26 + 3Ca(77)2/3 " Jq ^' 

That is 

(5-s)i/3>0(i)a(s)i/3. (89) 

Therefore, we may estimate the blowup rates at the origin of the density function p(s, x): 

lim p{s, 0){S - s)i/3 = lim -^{S - sf'^ > 0(1). (90) 

s^s- s-^s- a(s)'-''^ 

The proof is completed. ■ 

We notice that our blowup rate in the above theorem is different form the results of the 
Constantin and Escher's papers [71 and [1]. 
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